Universal Renyi mutual information in classical systems: the case of kagome ice 
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We study the Renyi mutual information of classical systems characterized by a transfer matrix. We first 
establish a general relationship between the Renyi mutual information of such classical mixtures of configuration 
states, and the Renyi entropy of a corresponding Rokhsar-Kivelson-like quantum superposition. We then focus 
on chiral and nonchiral kagome-ice systems, classical spin liquids on the kagome lattice. The chiral kagome ice 
is a critical phase realized in kagome layers of the pyrochlore spin ice in a magnetic field, while the nonchiral 
version has a finite correlation length and has been observed in hexagonal arrays of ferromagnetic nanoislands 
(artificial spin ice). Through a mapping of the chiral kagome ice to the quantum Liftshitz critical field theory, 
we predict a universal subleading term in the Renyi mutual information of this classical spin liquid. We verify 
our prediction with direct numerical transfer-matrix computations, and further demonstrate that the nonchiral 
kagome ice is a topologically trivial phase. 

PACS numbers: 



I. INTRODUCTION 

Universal subleading terms in the entanglement entropy en- 
code valuable information about quantum systems. In gapped 
topological phases, such as the Z2 spin-liquid ground state of 
the Kitaev toric-code model, 1 there is a subleading correc- 
tion to the area law, 2 known as the topological entanglement 
entropy, which encodes the total quantum dimension of the 
system. 3 4 This ground state of the toric-code model can be 
viewed as an equal- weight superposition of loops. A very sim- 
ilar universal correction can appear in critical wave functions 
such as the ground state of the Rokhsar-Kivelson (RK) 5 quan- 
tum dimer models; the so-called RK state is again an equal- 
weight superposition of dimer coverings. 6 9 In this latter case, 
the universal term has a different interpretation: it encodes 
information about the stiffness constant of the underlying crit- 
ical field theory, and therefore the exponent of critical corre- 
lations. 9 

The concept of topological order in classical systems has 
attracted considerable interest recently. 10 15 In particular, by 
using the notion of mutual information (a generalization of the 
entanglement entropy to mixed states), it was shown in Ref.fTOl 
that the topological order can survive decoherence and appear 
in classical systems. Moreover, it can be detected through 
subleading corrections to the mutual information. In particu- 
lar, it was argued in Ref. 10 that under fairly general condi- 
tions, the classical von Neumann mutual information is half 
of the corresponding quantum coherent case. Primarily focus- 
ing on critical systems, in this work we examine the universal 
terms in the Renyi mutual information of classical (thermal) 
mixture of configuration states, characterized by a transfer 
matrix, through a direct combinatorial approach. We obtain 
a general relationship between I c AB {n), the Renyi mutual infor- 
mation of the thermal mixture (where n is the Renyi index), 
and I q AB {n), of the corresponding quantum superposition: 
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magnetically disordered manifolds of spin configurations sat- 
isfying local constraints broadly referred to as ice rules. These 
systems can be thought of as the classical analogues of quan- 
tum spin liquids, whose entanglement properties (in particu- 
lar their topological entanglement entropy) have been the sub- 
ject of numerous recent studiesP^HHl D eS pke the prevalence 
of spin-ice systems in nature, and the connection to quantum 
spin liquids, the entanglement properties of spin-ice systems 
are largely unexplored. Here we focus on two spin-ice sys- 
tems on the kagome lattice, known respectively as chiral and 
nonchiral kagome ice. The local constraints, i.e., ice rules, 
characterizing these spin-ice manifolds allow for a transfer- 
matrix formulation of the problem, which in turn makes the 
general relationship ([T| applicable. By using a mapping of the 
chiral kagome ice to dimers on the honeycomb lattice, we re- 
late the universal correction to the Renyi mutual information 
of this classical spin liquid to universal terms in the Renyi en- 
tanglement entropy of the Liftshitz quantum critical point, and 
find 



Tab(») = 



1 



2(1 -n) 



In 



n + \ 



(2) 



Spin-ice 16 systems provide a canonical example of such 
classical thermal mixtures: they are massively degenerate 



where n is the Renyi index. Through a numerical transfer- 
matrix calculation, we also explicitly demonstrate that y c AB {n) 
vanishes for nonchiral kagome ice. 

Chiral kagome ice provides an exactly solvable experimen- 
tally relevant example, where subleading terms, analogous 
to the topological entanglement entropy, appear in the ab- 
sence of of quantum coherence or quantum fluctuations. The 
additional constraints in chiral kagome ice (with respect to 
the nonchiral manifold) result in critical spin correlations, as 
well as aforementioned subleading term y c AB {ri). In contrast, 
the nonchiral kagome ice, has a finite correlation length, ex- 
tremely weak dependence of the Renyi mutual information on 
the Renyi index (which likely vanishes in the thermodynamic 
limit), and a vanishing topological Renyi mutual information. 
Our general formulation in terms of a transfer matrix pro- 
vides a tool for studying the entanglement properties of more 
complex classical systems such as the three-dimensional py- 
rochlore spin-ice, as well as their RK-like quantum analogues. 

The outline of our paper is as follows. In Sec. Ill] we present 



the general setup of the problem, and derive Eq. ([TJ. We de- 
scribe in Sec.|III]the chiral and nonchiral kagome-ice systems, 
and discuss the relationship between the chiral manifold and 
the Liftshitz quantum critical point. We predict the universal 
subleading correction |2} using field-theory results on the Lift- 
shitz quantum critical point, and verify our prediction with di- 
rect numerical calculations. We also explicitly verify the van- 
ishing of the topological Renyi mutual information in nonchi- 
ral kagome ice. We close the paper in Sec. IV with a brief 
discussion. 
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FIG. 1: a) A convenient partitioning for studying universal entangle- 
ment is obtained by dividing an infinite cylinder of circumference I 
into two pieces, b) Folding the B subsystem onto the A half allows 
for the application of boundary conformal field theory. 



II. GENERAL SETUP 

The Renyi entropy Sx(n) of a subsystem X is a measure 
of entanglement with desired properties such as additivity and 
continuity: 
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Here the reduced density matrix px of the subsystem X is de- 
fined through tracing out the degrees of freedom in the rest 
of the system. The Renyi entropy depends on a parameter 
n known as the Renyi index, and potentially contains more 
information than the standard von Neumann entropy that is 
given by the n — > 1 limit of the Renyi entropy. 

For simply-connected subsystems, the knowledge of the 
Renyi entropies (for all indeces «) can yield the full spectrum 
of the the reduced density matrix, known as the entanglement 
spectrum, thus providin g a m ore complete characterization of 
entanglement properties p 3 l 24 l It has been shown, however, that 
in a wide class of gapped topological phases, the universal cor- 
rection to the Renyi entropy is independent of n. 25 For critical 
phases, on the other hand, the universal correction generally 
has a nontrivial dependence on n. 

The Renyi entropy is symmetric for a pure-state wave func- 
tion, i.e., Sa(h) = S b(k) for a system comprised of two sub- 
systems A and B. On the other hand, for generic mixed states 
we have Sa(h) + S b(«). A natural generalization of the Renyi 
entropy to mixed states, known as the mutual information 
7 AB (n), remedies this problem through an appropriate sym- 
metrization: 



Iab - j(Sa + S b - Saub), 



(4) 



where AU B represents the whole system, and we have sup- 
pressed the explicit dependence on the Renyi index n for 
brevity. Since, for a pure state, Saub(«) = for all n, the 
mutual information reduces to the Renyi entropy of entangle- 
ment in this case. The mutual information encodes the total 
amount of quantum and classical correlations between the two 
subsystems. 26 

For a system exhibiting area law with a universal order-one 
correction, we have^ 
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where c is a nonuniversal prefactor, £ is the d - 1 -dimensional 
"area" of the af-dimensional subsystem, y is the universal sub- 



leading term, and the dots represent other subleading contri- 
butions, which vanish in the limit of £ — > oo. In general, there 
could be logarithmic corrections to area law for subsystems 
with sharp corners or complex topology. 28 Here we only con- 
sider subsystems with vanishing logarithmic corrections. 

A convenient partitioning for extracting universal contribu- 
tions to the entanglement entropy [see Eq. |5])] is obtained 
by dividing an infinite cylinder of circumference £ in half™ 
The subsystems A and B then have simple topology and no 
sharp corners. We thus expect the area law to take the form of 
Eq. (BJ, with the subleading terms indicated by dots vanishing 
as e - ^ for gapped phases (where £ is the correlation length), 
and as a power-law of £ for critical systems. In practice, it is 
helpful to work with a finite cylinder of length 2/t (with free 
or fixed boundary conditions at the two endpoints), and take 
the thermodynamic limit h — > oo at the end of the calculation. 

Let us now consider a general set of (classical) configu- 
rations C satisfying some local constraints which allows for a 
transfer-matrix formulation. A Hilbert space corresponding to 
this manifold can be constructed as in the RK model: for each 
configuration we construct a state |C), and impose the orthog- 
onality condition (C\C) = Sec- If the degrees of freedom in 
a given configuration are classical Ising spins, then each state 
\C) is a direct product of single-spin up and down states | T) 
and | i). Assuming that there are N such configuration in a 
given (finite) system, an equal-weight quantum superposition 
and an equal-weight classical (thermal) mixture respectively 
have the following density matrices: 
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We start with the classical density matrix p c in Eq. |6|. Rep- 
resenting the degrees of freedom (say Ising spins) in subsys- 
tems A and B with J[ and S, we can write 
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where (I^IS)} constitutes an orthonormal basis. We further as- 
sume that the configuration states \C) belong to the set of basis 
states (which is certainly the case for Ising spins or other local 
degrees of freedom). We then find that the only contribution 
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FIG. 2: a) To calculate S C A , we need n replica states |C;), which have 
the same configurations in A subsystem (top panel). If we fold the 
B side onto the A side as in Fig. [T] (bottom panel), we get n + 1 
independent (except at the boundary) sets of configurations (or fields 
<pj in a field-theory representation) living on the half-cylinder, b) To 
calculate S^., we similarly identify the configurations on the B side, 
c) For the quantum case {S q A - S 9 B ), the configurations are identified 
with a different patter resulting in In independent fields in the bulk 
of the half-cylinder. 



to the sum above comes from |C,) = \3K-JBi) - ffli+ilii), which 
gives J{\ - 9\q_ — . . . J\. n . In other words, the n replica config- 
urations C, are "stitched" together as shown in Fig.[2k. 

Now if two configurations match in subsystem A, they must 
also match at the boundary between A and B. The sum in 
Eq. (|7]i above can then be written as a sum over the configura- 
tions d at the boundary as follows: 



nb ln (/^IX (d)i « 
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where NA(d) (NB(d)) denotes the number of configurations in 
A (B) for a fixed configuration d on the boundary (see Fig. [2^). 
The entropy S C B can be obtained simply by replacing A «-» 
B in the above expression (see Fig. |2t>). It is also easy to 



observe that S c , 



In N for all n. A similar treatment for 



the quantum density matrix p q in Eq. Q shows that the 2n 
configurations (C, and C\) are "stitched" together as shown in 
Fig.[2fc in the quantum case, and we can write 
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As the configurations C are characterized by local con- 
straints, the combinatorial problem of computing the entropies 
above can be formulated in terms of a transfer matrix, which 
has proved useful in the calculation of entanglement entropy 
in both static 7 and dynamical 29 cases. Let us break up the 
system into rows parallel to the boundary between the two 
subsystems, and indicate the configuration of the degrees of 
freedom on a row by a lower-case letter (the boundary con- 
figuration d indicates the configuration of a row at the inter- 
face of the two subsystems). In general, there are degrees of 
freedom between two consecutive rows (with configurations i 
and j), and the number of allowed configurations for these de- 
grees of freedom (for fixed i and j) can be denoted by 7y. We 
can now represent the row configurations by states and write 
Tij = (i\T\j), where T is the transfer matrix. 



Representing the states at the left- and right-hand side of 
the system, and the boundary between A and B (see Fig. [2^), 
respectively by \l), \r), and \d), we can then write 

N A {d) = ^</|r'V>, N B {d) = Yj(d\T h \r), (10) 

l r 

where we have assumed free boundary conditions at the left 
and right boundaries. Taking the thermodynamic limit, h — > 
oo, gives 

N A (d) = A h J](l\A)(A\d), N B (d) = A h J]<d|/l>(i|r>, 

/ r 

(ID 

where A is the largest eigenvalue of T, and \A) is the corre- 
sponding eigenvector. Similarly, we have N = Yjir(l\T 2h \r) - 
A 2h Z lr (l\A)(A\r). 

Note that for fixed (Dirichlet) boundary conditions (instead 
of free) we can use one particular \r) or \l) and omit the sum- 
mation over r and 1. Our final result will not change as the 
summation over the boundary conditions cancels out but we 
may need to choose a different \A): the largest-eigenvalue 
eigenstate of T, which has a nonvanishing overlap with both 
fixed \r) and \l). Also note that more generally, the density 
matrices may not represent equal-weight mixtures (or super- 
positions in the quantum case). As long as the configuration- 
dependent weights are local, they can be absorbed into the 
transfer matrix (in this case the number of configurations with 
given boundary conditions transforms into a partition function 
with those boundary conditions). 

Putting all these together, and assuming the overlaps be- 
tween \A) and the boundary states above are real, we then ob- 
tain 
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For this particular geometry, the expression for S C B is 
given by changing the dummy variable / to r (and yields 
an identical result). Similarly, we can write S C AUB = 

—^ In \A lh{ "- l) &i{A\l)f n ~ 2 \ which leads to the general ex- 
pression 



7 ^ (n) = TT^ ln Z w> " +1 



(13) 



where the Renyi index n is written explicitly. The entangle- 
ment entropy of the quantum case was studies in Ref. 7 using 
this transfer matrix method. For completeness, we restate that 
result in our notation: by inserting Eq. ( fTT) into Eq. (|9}, we 
immediately obtain 



7 !» = i^ ln [l> |d > 2 " 



(14) 



Notice that the row configurations d form an orthonormal 
basis and the normalization of the eigenvector \A) implies that 




FIG. 3: A kagome layer in the pyrochlore lattice. 



2 ( /(/l|<f ) 2 = 1 . In other words, p^ = (A\d) 2 is a probability, and 
the above expression is identical to that of Ref. 7 Comparing 
the two equations above, we finally arrive at the general rela- 
tionship I c AB (n) - \J^ B \r^-\ which is one of the main results 
of this paper [see Eq. {!}], and provides a nontrivial general- 
ization of the prediction of Ref. [10] 



III. KAGOME ICE AND QUANTUM LIFTSHITZ MODEL 

Highly frustrated magnets are canonical examples of clas- 
sical spin liquids with nontrivial correlations. In particular, 
the so-called spin-ice materials in which Ising-like moments 
reside on a three-dimensional network of corner-sharing tetra- 
hedra (the pyrochlore lattice shown in Fig. [3j exhibit a crit- 
ical dipolar-like spin correlations in an emergent Coulomb 
phase. 30 The origin of this critical phase can be traced to the 
nontrivial local constraints known as the two-in-two-out ice 
rules that govern the ordering of magnetic moments in indi- 
vidual tetrahedra. 16 

In this work, we apply the above analysis to two- 
dimensional analogues of the pyrochlore spin ice: (i) chiral 
and (ii) nonchiral kagome ice. The kagome lattice shown 
in Fig. Hk consists of corner-sharing triangles with two dif- 
ferent orientations. The projections of the easy axes on the 
kagome plane form 120-degree ordering as in Fig.|4^. In chi- 
ral kagome ice every up triangle has a two-in-one-out spin 
configuration with a net positive magnetic charge, while every 
down triangle has a one-in-two-out configuration with a net 
negative charge (see Fig. |4^); the entropy density of this ice 
phase is 0.108 k B per spin. 7132 Since the two types of triangles 
contain opposite magnetic charges, hence breaking the sublat- 
tice symmetry, this phase is referred to as the chiral kagome 
ice. The nonchiral kagome ice is a less constrained manifold: 
every triangle (whether up or down) is in either the two-in- 
one-out or the one-in-two-out state; the entropy density of this 
ice phase is 0.501 k B per spinPS. 

The nonchiral degenerate manifold minimizes the nearest- 
neighbor spin interaction energy. This ice phase is non- 
critical with short-ranged spin-spin correlations. The chiral 
kagome ice appears in the intermediate plateau regime of 
the pyrochlore spin-ice under a [111] field,^122l2il where the 
spins align themselves with the field on the triangular layers, 




FIG. 4: a) A chiral kagome ice configuration, b) The mapping to 
dimer coverings. The configurations of the bonds cut across by the 
sold black line represent a boundary state. 



and, consequently, form the degenerate chiral manifold on the 
kagome layers (see Fig. [3}. An artificial version of the kagome 
ice is also realized in arrays of single-domain ferro magnetic 
nanoisland arranged in a honeycomb networkP32Sl In artifi- 
cial kagome ice, the chiral ice phase is induced by further- 
neighbor spin-spin interactions 1221221 

Remarkably, spin-spin correlations in chiral kagome ice de- 
cay algebraically. 39 The critical correlations of this phase can 
be understood by mapping the chiral ice configurations to 
dimer coverings on the dual honeycomb lattice obtained by 
connecting the centers of the corner-sharing triangles of the 
kagome lattice 31 3239 (see Fig. Bp). Since every triangle has 
exactly one negative spin (lying on a bond of the honeycomb 
lattice), each site of the honeycomb lattice is visited by exactly 
one negative spin, which we identify with a dimer. The dimer- 
covering problem on honeycomb is exactly solvable using the 
well-known Pfaffian method for computing the partition func- 
tion. 40 Exact calculations showed a power-law dimer-dimer 
correlation decaying as l/r 2 F^-In the language of dimers, the 
nonchiral kagome ice then corresponds to dimer coverings on 
the honeycomb lattice with one type of defect: a site visited 
by exactly two dimers. 

In order to obtain the universal correction term to the Renyi 
mutual information for the chiral kagome ice [Eq. d2l], we 
first consider the Renyi entropy of the corresponding quan- 
tum RK state and then invoke the classical-quantum corre- 
spondence qTj). Universal entanglement entropy in RK wave 
functions of dimer coverings is a well-studied problem, 6 9 but 
the connection to spin-ice systems has not been exploited to 
study the the entanglement properties of these classical spin 
liquids. Let us briefly review the field-theoretical approach 
to the quantum RK case. The RK wave function \W) is an 
equal weight superposition of all dimer coverings on the hon- 
eycomb lattice, or equivalently an equal weight superposition 
of all chiral kagome-ice states. 

Dimer models on a bipartite lattice in turn have a height- 
model representation: scalar heights are defined on every 
hexagonal plaquette, when a small circular path around a 
site (which crosses three bonds) goes across an empty bond 
(dimer), we decrease (increase) the height by 1 (2). The 
path is directed clockwise for sites on one sublattice and 
counterclockwise for sites on the other. We then obtain a 
unique height field modulo the height of one reference pla- 
quette. A periodic structure is imposed on the heights to 
avoid overcounting configurations. 42 Upon coarse-graining, 



the height-model representation then leads to a critical field 
theory of compactified noninteracting bosons (coarse-grained 
height field) with effective actionPSBl 

S[0] = j- f c/ 2 x(d^) 2 , </>-(/> + 2nR, (15) 

which represents a universality class characterized by the dy- 
namical exponent z - 2, stiffness g, and the compactification 
radius R. By rescaling the bosonic fields <p, we observe that 
the only free parameter of the theory is R yfg, which can be 
computed by comparing the dimer correlation functions ob- 
tained from the above continuum quadratic filed theory with 
exact results based on the Pfaffian method.^ 

The projection of the RK wave function to a particular 
height field configuration and the quantum density matrix [see 
Eq. (|6]l] are given by 



<0|T> <* . 



-sw 



fft = i»p><n 



(16) 



The calculation of the entanglement entropy in the cylinder 
geometry then involves folding half of the cylinder (subsystem 
B) onto the other half. The identification of Fig. [2fc results in 
2n (n form each subsystem) independent (in the bulk) fields 
with say Dirichlet boundary condition on the left-hand side of 
the folded system (the far ends of both A and B subsystems), 
and a special boundary condition (p\ - <f>2 - ■ ■ ■ - <pi n on the 
right-hand side of the folded system (the boundary between 
the A and B subsystems). 

By explicitly constructing the boundary state correspond- 
ing to this boundary condition, it was shown in Ref. 9 (see 
also Refs. 17181 1 that the universal Renyi entropy should go as 

(N/2-1) 



\-n 
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where N is the number of inde- 



pendent fields, which in the quantum case equals N = 2«. 
Exactly the same argument also holds for the classical case 
except that, as demonstrated in Figs. [2^ and |2j), the num- 
ber of independent fields is AT = n + 1 (n from one sub- 
system and 1 from the other) with the boundary condition 
<p\ — <Pi - ' ' ■ = <f>n+\ on me right-hand side of the folded 
system. This is consistent with our Eq. ( fl3l > found through 
an explicit combinatorial approach. The actual stiffness of 
the dimer model corresponds to ~j2gR - 1 (see, e.g., Ref. 7), 
which immediately leads to our prediction [Eq. Q] of the uni- 
versal mutual information in chiral kagome ice. 

We also explicitly verified Eq. |2) with the transfer-matrix 
approach, i.e., by using Eq. ( p"3j ). We assume that the bound- 
ary between the A and B subsystems cuts across £ bonds of 
the honeycomb lattice (see Fig. Rb). The transfer matrix is 
then a 2 l x 2 f -dimensional matrix, which can be constructed 
using the local constraints (counting the number of dimer cov- 
erings by transfer matrices was pioneered by Lieb.^j). The 
transfer matrix on the honeycomb lattice has a block-diagonal 
structure: if the transfer matrix has a nonzero matrix element 
(i\T\j), the states |/) and \j) must have the same number of 
dimers. A mapping from dimers to fermions yields elegant 
analytical results for the transfer matrix both on the square 46 
and the honeycomb lattices. 7 

In the case of free boundary conditions at infinity, we need 
overlaps (A\i) for the eigenstate with the largest eigenvalue and 



all of the 2 f row states. As the transfer matrix is block diag- 
onal, we can consider only the states in the block with the 
largest eigenvalue. It was shown in Ref. 7 that the number of 
dimers in this block is £/3, Denoting the positions of dimers 
in state \i) as a, with < a\ < ai < ■ ■ • < cue/3 < £ — 1, the 
mapping to fermions allows us to compute |(/l|/)| 2 as a Van- 
dermonde determinant^ 



\(mr = 
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The above expression provides an efficient way to compute 
the Renyi entropy for large I (here we went up to t - 30). 
We have also checked the result by a brute-force method (nu- 
merical construction of the block-diagonal transfer matrix and 
direct exact diagonalization) for £ up to 18. The results are 
shown in Fig. Bh for n - 1 ... 4. As expected, for each n, the 
mutual information I c AB (n, €) satisfies the area law, i.e., scales 
as oc £ for large £. To extract the universal subleading term 
y, we fit this data to / = cl + y + d /■£, and obtained y as a 
function of n. The agreement with Eq. (|2} is excellent as seen 
in Fig. g). 

It is also illuminating to consider the case of nonchiral 
kagome ice, which has less constrained ice rules than the chi- 
ral manifold. The analogous quantum wave function then has 
a finite correlation length. The subleading correction to the 
area law in this case can then be interpreted as topological 
entanglement entropy. We find through numerical transfer- 
matrix calculations that the mutual information satisfies the 
area law, and a fit to cl + y gives y of order 1CT 4 , i.e., the 
nonchiral kagome ice (and the corresponding quantum RK 
wave function) is topologically trivial. The transfer matrix 
does not have a block-diagonal structure in this case, and we 
do not have an exact analytical solution for the overlaps (A\i) 
so we are limited to smaller systems in the numerics. How- 
ever, we do not need to go to larger systems because the finite- 
size corrections are exponentially small. The universal terms 
can then be easily extracted from a simple linear fit for very 
small system sizes. 

An interesting question is the behavior of y when we in- 
terpolate between these two manifolds. If we consider the 
pyrochlore spin ice in a magnetic field H in the [111] direc- 
tion, we have a density matrix that only depends on the ratio 
H/T at temperature T. In terms of the dimer model, there is 
an energy cost proportional to H for any defect (two dimers 
touching on one site). At zero temperature (and H > 0), 
H/T diverges and no defects are allowed, which results in 
the critical chiral kagome-ice manifold. On the other hand, 
the nonchiral ice corresponds to H/T = 0, where there is no 
cost for such defects. At low temperatures, any infinitesimal 
finite temperature destroys the critical phase, resulting in a 
finite density of defects and a correlation length that is expo- 
nentially large in H/T? 1 This instability follows from the fact 
that in the height representation, a defect is like a vortex (the 
compactified height field winds once when going around such 
defect) with an operator proportional to cos ( y26), where 6 
is the dual field to the height field (p, which is a relevant 
(in renormalization-group sense) perturbation. 47 The fate of y 
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FIG. 5: a) Numerically computed mutual information using Eqs. l(T7} 
and ([13} exhibits the area law. b) Fitting the data in panel (a) gives 
the universal piece y (blue stars), which are in perfect agreement 
with the analytical prediction Eq. ([2} (red line), c) For the nonchiral 
kagome ice, the mutual information has a very weak dependence on 
the Renyi index n, and the linear fit to data gives y c AB a 0. 



[see Eq. Q] at small finite temperatures (large finite H/T) re- 
mains an open question as the explicit computation of y is dif- 
ficult for large (but nondivergent) H/T, due to the absence of 
an exact solution (in the presence of defects) and large finite- 
size effects (an exponentially large correlation length). 

We note, however, that this finite density of defects does 
not break the Zi chirality order parameter (similar to domain 
walls in an Ising ferromagnet), and this low-temperature non- 
critical phase is distinct from the nonchiral ice, where the Z2 
symmetry is unbroken. Despite a finite correlation length, it 
seems likely that Eq. (|2} should remain unchanged as we in- 



crease temperature. At some critical temperature, we expect 
an Ising-like phase transition to the nonchiral phase, where 
y should jump to zero. In fact for small H/T, where we 
have a small correlation length, we have performed numeri- 
cal transfer-matrix calculations, and verified the presence of a 
y - plateau as we move away from the nonchiral kagome 
ice corresponding to H/T = 0. As the correlation length keeps 
increasing with H/T, however, the calculation can not be ex- 
tended to the vicinity of the the chiral phase. 



IV. DISCUSSION 

Motivated by the predictions of Ref. 10 on the appear- 
ance of topological entanglement entropy in classical systems, 
which are obtained through the decoherence of topologically 
ordered quantum wave functions, we obtained a generalized 
relationship [Eq. ([TJ] between the Renyi mutual information 
of generic classical systems (which can be described by a 
transfer matrix), and their quantum-coherent RK-like coun- 
terparts. The classical systems above can be similarly viewed 
as decohered versions of these RK quantum wave functions. 

We examined the universal entanglement properties of two 
spin-ice manifolds (chiral and nonchiral) on the kagome lat- 
tice, and found a nonvanishing universal subleading term in 
the mutual information of the chiral state (a highly constrained 
classical spin liquid with critical correlations). The experi- 
mentally relevant chiral kagome ice, which appears in kagome 
layers of pyrochlore spin-ice in the presence of a [11 In- 
direction magnetic field, readily maps onto a dimer model on 
the honeycomb lattice. The RK quantum dimer model, and 
its well-known corresponding quantum Liftshitz universality 
class, thus describe the quantum-coherent counterpart of the 
chiral kagome ice. This mapping, and the general relation- 
ship ([TJ thus leads to our prediction of the presence of a uni- 
versal subleading term in chiral kagome ice, which bears a 
striking resemblance to the topological entanglement entropy 
in gapped topologically ordered phases. 

Unlike topological phases, where all Renyi entropies are the 
same, the universal mutual information in this critical classi- 
cal spin liquid, has a nontrivial dependence on the Renyi in- 
dex. Moreover, we verified that for the short-range correlated 
nonchiral kagome ice, the dependence of the mutual infor- 
mation on the Renyi index is extremely weak (likely only in 
the subleading terms), and the topological mutual information 
vanishes. 
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